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Me  have  made  significant  progress  this  year  in  two  new  program  endeavors 
which  will  impact  all  areas  of  our  ONR  Power  Programs  research  on  energetic  com¬ 
pounds. 

I.  We  have  so  far  completed  the  analysis,  coding  and  testing  of  the  inte¬ 
gral  and  SCF  portions  of  the  crystal  orbital  part  of  the  program  and  the  analy¬ 
sis  and  some  coding  for  the  polymer  orbital  part.  We  meshed  into  this  new  POLY- 
CRYST  program,  as  options,  all  of  our  desirable  computational  strategies  for  ab- 
initio  calculations  on  large  molecules. 

II.  Ab-initio  coupled-cluster  program  including  multi  reference  determinants 
and  a  new  original  repetitive  method,  based  either  on  localized  orbitals  or  con¬ 
ventional  delocalized  canonical  orbitals. 

Our  strategy  in  this  part  of  the  research  has  been  to  mesh  together  with 
configuration  interaction  the  desirable  features  of  coupled-cluster  theory.  We 
have  carried  out  preliminary  calculations  on  RDX. 

III.  We  have  continued  optimizing  geometries  and  generating  electrostatic 
molecular  potential  contour  maps  (EMPC)  maps  for  new  monomers  of  interest  ex¬ 
perimentally. 

This  past  year  we  have  also  been  working  to  develop  methodology  for  more 
quantitative  predictions  of  relative  polymerization  ratios.  Thus,  we  have  de¬ 
veloped  and  implemented  methods  to  calculate  the  volume,  shape  and  accessibil¬ 
ity  of  isopotential  EMPC  contours,  and  are  developing  a  strategy  for  comparing 
el ectrophili cities  for  the  propagation  step. 

IV.  Using  our  CRYSTAL  program,  electrostatic  and  electric  field  terms  based 
on  the  molecule  RDX,  and  atom  class-atom  class  potentials  from  CH3NO2,  and a 
series  of  diatomic  dimers,  we  have  predicted  the  density  of  RDX  to  within  0.1%. 
This  method  then  holds  promise  for  predicting  unit  cell  parameters  of  large  ex¬ 
plosives. 

V.  We  have  continued  to  carry  out  ab-initio  calculations  on  polynitropoly- 
hedranes  and  polynitropolyazapolyhedranes  to  generate  electrostatic  molecular 
potential  contour  (EMPC)  maps  around  these  polynitropolyhedranes  and  polynitro- 
jolyazapolyhedranes.  These  EMPC  maps  appear  to  indicate  that  the  character  of 
the  nitro  groups  may  change  as  a  function  of  the  degree  of  nitration,  and  they 
also  may  prove  useful  in  understanding  how  these  molecules  will  pack  in  cry¬ 
stals. 

VI.  Ab-initio  configuration  interaction  calculations  using  multi  reference 
determinants  have  been  carried  out  for  much  of  the  CH3-NO2  decomposition  path¬ 
way  of  nitromethane.  For  each  point,  in  the  order  of  800,000  configurations 
were  generated  and  sifted  using  perturbation  theory  for  the  most  important. 

VII.  Preliminary  Cl  calculations  were  carried  out  on  RDX  at  equilibrium  con¬ 
sidering  the  molecular  orbitals  centered  primarily  on  a  r 

c>-«i 

portion  of  the  skeleton  as  a  prototype  for  the  N-NO2  decomposition.  Even  at 
equilibrium,  a  multi  reference  determinant  character  on  the  NO2  group  was 
found— similar  to  that  in  CH3-NO2'. 

VIII.  We  have  gotten  access  to  CYBER  205  and  CRAY  time.  We  are 
the  process  of  converting  our  programs  to  both  of  these  machines. 
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We  have  made  significant  progress  this  year  in  two  new  program 
endeavors  which  will  impact  all  areas  of  our  ONR  Power  Programs  research 
on  energetic  compounds. 

I.  Ab-initio  crystal  orbitals  and  polymer  orbitals 

This  program  will  enable  one  to  answer  such  questions  as: 

How  do  the  occupied  and  virtual  molecular  orbital  energy  levels  and  the 
bond  strengths  of  a  molecule  change  from  the  isolated  molecule  when  the 
molecule  is  in  a  crystal? 

What  is  the  bonding  between  molecules  in  a  crystal? 

These  both  affect  properties  important  in  explosives. 

How  does  the  bonding  change  within  an  energetic  pendant  group  and  be¬ 
tween  the  pendant  group  and  the  backbone  between  a  monomer  and  a  poly¬ 
mer? 

We  have  so  far  completed  the  analysis,  coding  and  testing  of  the  integral 
and  SCF  portions  of  the  crystal  orbital  part  of  the  program  and  the 
analysis  and  some  of  the  coding  for  the  polymer  orbital  part  of  the 
program.  We  meshed  into  this  new  POLY-CRYST  program,  as  options,  all 
of  our  desirable  computational  stategies  for  ab-initio  calculations  on 
large  molecules. 

II.  Ab-initio  coupled  cluster  program  including  multi  reference  deter¬ 
minants  and  a  new  original  repetitive  method,  based  either  on  localized 
orbitals  or  conventional  delocalized  canonical  orbitals. 

Our  large-scale  ab-initio  multi  reference  determinants  configuration 
interaction  calculations  on  the  decomposition  of  nitromethane  CH.-N0? 
decomposition  pathway  indicated  that  not  only  were  there  the  two  deter¬ 
minants  found  at  the  equil ibrium  geometry  (located  on  the  NCL  group), 
but  also  that  as  the  C-N  bond  was  broken  additional  determinants  became 
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important. 


Our  strategy  in  this  part  of  the  research  has  been  to  mesh  together 
with  configuration  interaction  the  desirable  features  of  coupled  cluster 
theory. 

Thus,  for  calculations  of  decomposition  pathways  of  energetic  mol¬ 
ecules  with  NO-  groups,  multi  reference  determinants  will  be  necessary 
within  any  many-body  perturbation  theory  or  coupled  cluster  theory. 

We  have  completed  the  analysis,  coding  and  testing  of  a  coupled 
cluster  program  based  either  on  localized  orbitals  or  delocalized  or¬ 
bitals  and  have  carried  out  preliminary  calculations  on  RDX. 

We  have  so  far  completed  the  analysis,  coding  and  testing  for  the 
linear  coupled  cluster  terms  based  on  multi  reference  determinants  and 
have  carried  out  preliminary  calculations  on  RDX.  We  are  now  working 
on  the  implementation  of  the  higher  order  terms.  However,  since  we 
can  include  any  number  of  multi  reference  determinants  in  our  starting 
wave  function,  this  should  ameliorate  the  necessity  for  many  or  most 
of  the  higher  order  terms.  We  will  be  testing  this  hypothesis. 


Another  strategy  is  to  apply  the  T,  and  T-  operators  repetitively, 
each  time  on  the  improved  wave  function  from  the  previous  cycle.  This  I 

will  generate  the  higher  order  terms  but  has  the  advantage  of  keeping 

down  the  size  of  the  equations  to  be  solved.  We  have  completed  the 

analysis,  coding  and  testing  of  the  strategy  and  have  carried  out  pre¬ 
liminary  calculations  on  RDX. 

We  have  also  continued  our  research  and  made  gratifying  progress  • 

in 

III.  Ab-initio  quantum  chemical  calculations  and  electrostatic  molecular  « 

potential  contour  maps  for  cationic  polymerization  i 

We  have  continued  optimizing  geometries  and  generating  electro-  j 

static  molecular  potential  contour  (EMPC)  maps  for  new  monomers  of 

interest  experimentally.  Our  previous  results  had  illustrated  vividly  ] 

that  these  EMPC  maps  were  capable  of  predicting  reliably  the  basicity  : 

of  even  hypothetical  monomers  prior  to  their  synthesis.  The  initiation  1 

step  in  the  cationic  polymerization  of  cyclic  ethers  is  governed  by  the  j 

basicity  of  these  monomers  and  our  EMPC  maps  predicted  reliably  the  j 

rank  order  of  polymerization  and  even-  optimal  copolymer  candidates. 


This  past  year  we  have  also  been  working  to  develop  methodology 
for  more  quantitative  predictions  of  relative  polymerization  ratios. 
This  is  complicated  by  the  fact  that  there  are  two  steps  in  the  poly¬ 
merization:  initiation— proportional  to  basicity  and  propagation- 

considered  by  experimental  polymer  chemists  to  be  influenced  by  the 
electrophilicity  of  the  propagating  ring  opened  carbocation.  We  had 
already  shown  that  these  two  steps  are  influenced  in  opposite  and  com¬ 
peting  directions  by  the  character  of  the  substituent  groups.  We  had 
also  shown  that  initiation  was  a  dominant  step.  Thus  we  have  developed 


I 
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and  implemented  methods  to  calculate  the  volume,  shape  and  accessibility 
of  isopotential  EMPC  contours.  We  are  now  carrying  out  such  calculations 
to  compare  against  experimental  data.  We  are  also  developing  a  strategy 
for  comparing  electrophil icities  for  the  propagation  step. 

IV.  Prediction  of  crystal  density 

We  have  been  quite  successful  in  this  endeavor  during  the  past  year. 

Using  our  CRYSTAL  program,  we  predicted  the  unit  cell  volume  of 
CH3N02  to  within  4.27%  and  hence  the  density  to  within  -4.2 7%  using  ab- 
inTtio  atom  class  -  atom  class  potential  functions  from  energy  parti¬ 
tioned  ab-initio  calculations  on  CkLNCL  plus  short  range  atom  -  atom 
terms  from  a  series  of  diatomic  dimers; 

Certain  of  the  terms  such  as  themultipole  contributions  to  the 
electrostatic  term  and  the  electric  field  for  the  polarization  term  are 
recalculated  explicitly  each  time  for  each  change  in  unit  cell 
parameter  from  the  wave  function  of  the  molecule. 

Even  more  gratifying  is  that  we  predicted  the  unit  cell  dimensions 
of  RDX  to  within  -1.0  (a),  1.5  (b),  -0.6  (c)%  and  thus  the  density  to 
within  0.1%,  employing  the  same  strategy  but  using  the  ab-initio  atom 
class-atom  class  potentials  from  CH3NO2  and  the  diatomic  dimers. 

This  then  holds  promise  for  predicting  unit  cell  parameters  of  large 
explosives.  Such  calculations  are  very  sensitive  to  the  goodness 
of  the  potentials. 

V.  Ab-initio  calculations  on  polynitropolyhedranes  and  polynitropolyaza- 

polyhedranes 

We  have  continued  to  carry  out  ab-initio  calculations  on 
polynitropolyhedranes  and  polynitropolyazapolyhedranes.  From  these 
wave  functions  we  have  been  generating  electrostatic  molecular  potential 
contour  (EMPC)  maps  around  these  polynitropolyhedranes  and  polynitro¬ 
polyazapolyhedranes.  These  EMPC  maps  appear  to  indicate  that  the 
character  of  the  nitro  groups  may  change  as  a  function  of  the  degree  of 
nitration.  Also  these  EMPC  maps  may  prove  useful  in  understanding  how 
these  molecules  will  pack  in  crystals. 

Molecules  which  have  a  number  of  nitro  groups  around  their  peri¬ 
pheries  could  scarcely  be  expected  to  approach  each  other  closer  than 
a  zero  isopotential  contour  if  the  signs  of  the  EMPC  contours  of  the 
adjoining  molecules  are  the  same.  On  the  other  hand,  if  the  molecules 
are  surrounded  by  molecules  with  both  positive  and  negative  EMPC  contours, 
the  molecules  may  try  to  orient  themselves  in  a  crystal  having  regions 
of  positive  EMPC  contours  around  one  molecule  adjoining  regions  of 
negative  EMPC  contours  around  the  other  molecule. 
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VI.  Ab-initio  configuration  interaction  calculations  on  the  CH3-NO2 
decomposition  pathway  of  nitromethane 

Ab-initio  configuration  interaction  calculations  using  multi  reference 
determinants  have  been  carried  out  for  much  of  the  CH3-NO2  decomposition 
pathway  of  nitromethane.  For  each  point  in  the  order  of  800,000 
configurations  were  generated.  Using  a  perturbation  procedure  the  energy 
contributipn  of  each  configuration  was  calculated.  All  configurations 
contributing  an  energy  higher  than  a  chosen  threshold  were  included  in 
the  Cl  wave  function  to  be  solved.  Then  an  extrapolation  was  performed 
to  take  into  account  the  energy  contributions  of  the  configurations  not 
included  explicitly.  This  was  then  followed  by  a  Davidson  type 
correction  for  size-consistency. 

VII.  Ab-initio  configuration  interaction  calculations  on  RDX 

Preliminary  Cl  calculations  were  carried  out  on  RDX  at  equilibrium 
considering  the  molecular  orbitals  centered  primarily  on  a 


portion  of  the  skeleton  as  a  prototype  for  the  N-NO2  decomposition. 
Even  at  equilibrium  a  multi  reference  determinant  character  on  the  N02 
group  was  found  -  similar  to  that  in  CH3-NO2.  Thus,  it  appears  to 
be  a  general  characteristic  of  all  aliphatic  nitroexplosives. 


Several  major  invited  papers  have  been  presented  on  our  polymer 
research  at  the  1983  Gordon  Conference  on  Polymers  and  at  the  1983 
NATO  Advanced  Study  Institute  on  Quantum  Theory  of  Polymers:  Solid 
State  Approaches.  Our  plenary  lecture  at  the  NATO  Advanced  Study  Institute 
also  included  results  using  our  CRYSTAL  program  on  prediction  of 
crystal  densities.  An  invited  paper  on  the  crystal  results  was  also 
presented  at  the  Sanibel  International  Symposium  on  Molecular  and  Solid 
State  Theory  and  Computational  Quantum  Chemistry. 

Papers  on  the  polymer,  crystal,  polynitropolyhedrane  and  nitromethane  Cl  results 
were  also  presented  at  the  National  American  Chemical  Society  and 
American  Physical  Society  Meetings. 

Papers  on  several  of  the  above  topics  are  already  accepted  for 
publication  and  in  press. 


VIII.  Conversion  of  our  programs  to  vector  supercomputers,  CYBER  205 
and  CRAY 

The  quantum  chemical  problems  we  are  investigating  involve  large 
scale  calculations  on  large  systems,  both  with  many  atoms  and/or  with 
many  configurations  or  simultaneous  equations  to  be  solved.  We  have 
gotten  access  to  CYBER  205  time  and  to  CRAY  time.  The  CYBER  205  has 
two  million  64  bit  words  in  real  memory  plus  virtual  memory.  The 
CRAY  has  one  million  words  of  real  memory.  We  are  in  the  process  of 
converting  our  programs  to  both  of  these  machines. 


QUANTUM  CHEMICAL  INVESTIGATIONS  OF  THE  MECHANISM  OF 
CATIONIC  POLYMERIZATION  AND  THEORETICAL  PREDICTION  OF  CRYSTAL 
DENSITIES  AND  DECOMPOSITION  PATHWAYS  OF  ENERGETIC  MOLECULES 

Joyce  J.  Kaufman,  Principal  Investigator 


I.  Ab-initio  crystal  orbitals  and  polymer  orbitals 

II.  Ab-initio  coupled-cluster  program 

III.  Ab-initio  quantum  chemical  calculations  and  electrostatic  molecular 
potential  contour  maps  for  cationic  polymerization 

IV.  Prediction  of  crystal  density 

V.  Ab-initio  calculations  of  polynitropolyhedranes  and  polynitropoly- 
azapolyhedranes 

VI.  Ab-initio  configuration  interaction  calculations  on  the  CH,-N02 
decomposition  pathway  of  nitromethane  J  c 

VII.  Ab-initio  configuration  interaction  calculations  on  RDX 

VIII.  Conversion  of  our  programs  to  vector  supercomputers,  CYBER  205  and 
CRAY 


We  have  made  significant  progress  this  year  in  two  new  program 
endeavors  which  will  impact  all  areas  of  our  ONR  Power  Programs  research 
on  energetic  compounds. 

I.  Ab-initio  crystal  orbitals  and  polymer  orbitals 

We  have  so  far  completed  the  analysis,  coding  and  testing  of  the 
integral  and  SCF  portions  of  the  crystal  orbital  part  of  the  program 
and  the  analysis  and  some  of  the  coding  for  the  polymer  orbital  part 
of  the  program.  We  meshed  into  this  new  POLY-CRYST  program,  as  options, 
all  of  our  desirable  computational  strategies  for  ab-initio  calculations 
on  large  molecules. 

A.  Methodology 

1.  Major  Features 

This  technique  permits  one  to  calculate  ab-initio  quantum  chemical 
crystal  orbitals  and  polymer  orbitals  making  use  of  the  translational 
symmetry  in  a  crystal  and  translational  and/or  translational -rotational 
symmetry  in  a  polymer. 
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Main  features  of  the  ab-initio  SCF  portion  of  the  crystal  orbital 
part  of  the  program  method: 

Fully  Ab-inito 

Full  use  of  translational  symmetry 
Analogous  to  molecular  Hartree-Fock-Roothaan  Method 
No  theoretical  problems  in  approaching  Hartree-Fock  lim¬ 
it  (All  numerical  difficulties  can  have  arbitrarily 
strict  convergence  criteria) 

General  in  1,  2  or  3  dimensions 

Differences  from  the  ab-initio  SCF  molecular  method: 

Periodic  Gaussians  replace  Gaussians  in  basis  set 
Additional  quantum  number  present  in  orbital  description 
Gross  atomic  orbital  populations  and  total  overlap  popu¬ 
lations  in  molecule  reflect  influence  of  surround¬ 
ing  molecules. 

2.  Problem  Definition 

Specify  the  location  of  the  molecules  within  one  unit  cell.  Spe¬ 
cify  a  set  of  lattice  vectors  {R}  which  describe  the  translation  of 
the  primary  unit  cell  to  all  equivalent  positions.  Include  enough  unit 
cells  in  the  integrals  until  the  interactions  are  negligible.  Calcula¬ 
tion  of  the  wavefunction  within  one  unit  cell  is  now  sufficient  since 
all  unit  cells  are  equivalent.  The  cell-cell  interactions  are  already 
included  in  the  integrals. 

a.  Introduction  of  k  into  the  Problem 

As  an  electron  moves  in  the  periodic  potential  of  the  crystal,  the 
difference  between  unit  cells  is  only  due  to  the  phase  difference. 

This  can  be  introduced  into  our  formalism  along  with  the  periodicity  by 
introducing  phase  factors  of  e-ik»R  into  our  functions.  It  is  not  con¬ 
venient  to  deal  directly  with  a  continuous  variable  in  our  calculations 
In  practice  we  evaluate  the  crystal  properties  on  a  mesh  of  discrete  £ 
points  which  can  be  as  fine  as  necessary  to  display  all  the  features  of 
the  band  structure. 


b.  New  Quantum  Number  $ 


It  is  well  know  that  bringing  identical  atoms  together  causes  a 
splitting  of  the  degeneracy  of  the  orbital  levels  in  the  formation  of 
molecular  orbitals.  Similarly,  assembling  identical  molecules  causes 
a  "smearing  out"  of  the  molecular  orbital  levels  into  energy  bands,  as 
the  energies  of  the  crystal  orbitals  are  usually  known.  The  position 
of  an  electron  within  the  band  is  described  by  the  continuous  variable 
k,  representing  the  "crystal  momentum"  or  "wavenumber"  of  the  electron. 
Because  of  the  translational  symmetry,  £  can  be  restricted  without  loss 
of  generality  to  the  first  Brillouin  Zone  (first  unit  cell  in  momentum 
space. ) 

The  form  of  our  periodic  basis  function  is  now: 

Ie‘1k‘R  x  (nR) 

■Jk. 

R 

where  R  runs  in  principle  over  all  lattice  vectors  of  the  crystal. 

Only  cells  with  interactions  above  our  threshold  must  be  included. 

For  xp ( r )  we  can  choose  any  convenient  functions,  such  as 
Gaussian-  or  Slater-type  orbitals.  We  choose  a  sum  of  primitive 
cartesian  Gaussians. 

The  form  of  the  answer  is 


♦nk(r) 


1 


N  B 

vff  ^  l 
R  w 


-ik-R 


'nu 


kv 


xp(r+R) 


where  N  is  the  number  of  cells  and  B  is  the  size  of  the  basis  set  in 
the  primary  unit  cell.  The  £  dependence  is  indicated  by  subscripts  on 
4>  and  c.  The  variation  within  £  is  smooth.  This  band  holds  2N  elec¬ 
trons,  so  the  normalization  of  /Jr  gives  the  correct  number  per  unit 
cell.  ^ 
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c.  Hartree-Fock  crystal  orbital  equations  and 
matrix  elements 

The  Hartree-Fock  crystal  orbital  equations  are 


ffnlt  tnt  =  EnTt 


k  unlt 


in  matrix  form.  The  overlap  matrix  elements  are 


"  -ilt-K  » 

V*  - 1 e  V1 


and  the  Hamiltonian  matrix  elements  are 


<  %<?> 
ft 


-v 
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xjr+ft): 


cell  N  ^  -ZT 

*  H  <  xM(?)  I  -Hr: 

i  #'  lei  -  rl 


xv(M)> 


N  N 


*1  l  l  i  /*’  I 

*  rj  y6  Brillouin  n 
Rl  R2  Zone 


.  c  x 


c  nyft'  n6ft 


[2<xy(?'-S2)  x„(?)  |  -f— 
'  M  I r ' -r  I 


Xy^J  x6(r'+fti): 


<xv(r'+^2)  x„(r)  |  I  X,s(r+^1)  X  ,(r'4)>]} 


I  r'  -r 


3.  Logic  of  Crystal  Orbital  Program 

Step  1.  Generation  of  lattice  vectors  from  cell  parameters  and  cell 
selection  criteria. 

Step  2.  Calculation  of  real-space  one-electron  integrals  for  all  cases 
where  one  function  is  in  the  primary  unit  cell. 

Step  3.  Calculation  of  real-space  two-electron  integrals  for  all  cases 
where  one  function  is  in  the  primary  unit  cell. 

Step  4.  For  each  k-point,  obtain  the  integrals  for  the  periodic  basis 
functions.  The  pieces  from  different  cells  must  be  collated 
and  phase  factors  introduced. 

Step  5.  Input  initial  coefficients  and  k-occupation.  / 

Steps  6-9  are  SCF  cycle 

Step  6.  Form  density  matrix  by  3-dimensional  numerical  integration 
over  £. 

Step  7.  Form  Fock  matrix  for  each  k. 

Step  8.  Diagonalize  each  Fock  matrix. 

Step  9.  Obtain  new  coefficients,  test  for  convergence  and  do  bookkeeping. 


4.  Remaining  Program  Development  and  Implementation 

a.  Population  Analysis 

We  have  done  the  analyses  and  are  now  implementing  the  population 
analysis  portion.  The  results  of  the  population  analyses  gives  the 
gross  atomic  charges  on  the  atoms  and  the  total  overlap  popualtional 
TOP'S  between  atoms  (both  intra  and  intermolecular).  These  TOP'S  in¬ 
dicate  whether  the  intramolecular  bonding  is  strengthened  or  weakened 
when  a  molecule  is  placed  in  a  crystal  or  a  polymer.  These  TOP'S  also 
indicate  the  strength  of  intermolecular  bonding. 

b.  Properties  Package 

The  analyses  has  been  done  for  a  properties  package  to  calculate 
atomic  and  molecular  multi  pole  moments.  These  are  important  in  inter¬ 
molecular  interactions.  This  section  will  be  implemented. 

c.  Inclusion  of  Correlation  Energy 

We  are  now  giving  careful  consideration  into  how  best  to  include 
correlation  effects.  At  the  NATO  Advanced  Study  Institute  on  Quantum 
Theory  of  Polymers:  Solid  State  Approaches  last  sunnier,  it  was  suggest¬ 
ed  that  using  many-body  perturbation  theory  or  coupled  cluster  theory 
based  on  localized  orbitals  might  be  the  most  tractable,  fruitful  ap¬ 
proach. 

We  had  already  coded  and  tested  for  molecules  the  coupled  cluster 
approach  based  either  on  localized  orbitals  or  on  delocalized  orbitals 
in  connection  with  our  research  on  the  decomposition  of  the  nitroex- 
plosives.  (See  next  section  II).  We  are  now  looking  into  how  best  to 
mesh  this  approach  with  the  crystal  orbital  /  polymer  orbital  program. 


Major  credit  for  this  effort  is  due  to  Dr.  John  M.  Blaisdell  with  the 
collaboration  of  Dr.  P.C.  Hariharan. 


B.  Types  of  Pertinent  Questions  Which  Can  Be  Answered  With 

This  Method 

This  program  will  enable  one  to  answer  such  questions  as: 

1.  Crystals 

How  do  the  occupied  and  virtual  molecular  orbital  energy  levels 
and  the  bond  strengths  of  a  molecule  change  from  the  isolated  molecule 
when  the  molecule  is  in  a  crystal? 

What  is  the  bonding  between  molecules  in  a  crystal? 

2.  Polymers 

How  does  the  bonding  change  within  an  energetic  pendant  group 
and  between  the  pendant  group  and  the  backbone  between  a  monomer  and 
a  polymer? 


12 


II.  Ab-initio  coupled  cluster  program,  including  multi  reference  determinants 

and  a  new  original  repetitive  method,  based  either  on  localized  orbitals 

or  conventional  delocalized  canonical  orbitals. 

A.  Introductory  discussion 

Our  strategy  in  this  part  of  the  research  has  been  to  mesh  together  with 
configuration  interaction  the  desirable  features  of  coupled-cluster  theory. 

Our  large-scale  ab-initio  multi  reference  determinant  configuration 
interaction  calculations  on  the  CH3-NO2  decomposition  pathway  of  nitro- 
methane  indicated  that,  not  only  were  there  the  two  determinants  found  at 
the  equilibrium  geometry  (located  on  the  NO2  group),  but  also 
C-N  bond  was  broken,  additional  determinants  became  important. 

Thus,  for  calculations  of  decomposition  pathways  of  energetic  molecules 
with  NO2  groups,  multi  reference  determinants  will  be  necessary  within  any 
many-body  perturbation  theory  or  coupled-cluster  theory. 

We  have  completed  the  analysis,  coding  and  testing  of  coupled-cluster 
theory  based  on  multi  reference  determinants. 

We  have  completed  the  analysis,  coding  and  testing  of  a  coupled-cluster 
program  based  either  on  localized  orbitals  or  delocalized  orbitals  and 
have  carried  out  preliminary  calculations  on  RDX. 

Another  strategy  is  to  apply  the  T-|  and  T2  operators  repetitively,  each 
time  on  the  improved  wave  function  from  the  previous  cycle.  This  will 
generate  the  higher  order  terms  but  has  the  advantage  of  keeping  down  the 
size  of  the  equations  to  be  solved.  We  have  completed  the  analysis,  coding 
and  testing  of  the  strategy  and  have  carried  out  preliminary  calculations 
on  RDX. 


B.  Background  and  methodology 
1.  Coupled-cluster  method 

The  coupled-cluster  method  is  an  attempt  to  introduce  interactions 
among  electrons  within  clusters  as  well  as  coupling  amoung  these  clusters 
of  electrons  and  to  permit  the  wave  functions  to  contain  all  possible 
disjoint  clusters.  Coupled-cluster  theory  had  previously  been  shown  by 
others  to  be  a  useful  approach  to  the  correlation  problems  in  small  and 
medium  size  atomic  and  molecular  systems.  The  coupled-cluster  method 
truncated  to  first  and  second  order  had  been  used  by  others  in  closed  and 
open  shell  cases. 

The  question  arises,  "How  important  are  higher  orders,  specifically 
third  and  fourth  order,  in  the  coupled-cluster  method?"  Previous 
studies  of  closed-shell  systems  had  demonstrated  that  fourth  order  terms 
are  often  chemically  significant. 

First,  the  full  set  of  diagrams  for  closed  and  open  shell  systems, 
including  all  single  and  double  excitations,  was  derived  originally 
by  Saut4  in  France  (for  his  doctoral  thesis  with  Professor  Laforgue). 


Then  in  our  group  at  the  Johns  Hopkins  University  these  diagrams  were 
systematically  tabulated  by  Professor  Laforgue.  Now,  in  our  group 
Professor  Laforgue  and  Or.  Roszak  have  expanded  coupled  cluster  theory 
through  third  and  fourth  order,  including  sets  of  equations,  complete 
sets  of  diagrams  and  the  table  of  matrix  elements.  This  latter  has  been 
implemented  computationally  at  the  Johns  Hopkins  University  by  Dr.  Roszak. 


In  the  cluster-type  expansion  the  exact  wave  function  is  expressed 
in  the  following  form: 

|  y>  =  eT|  <j>  >  (II-l ) 

where  |  v>  represents  the  wave  function  in  the  one-particle  model; 


T  is  the  wave  operator  and  e^  is  the  correlation  operator, 
generates  one-,  two-electron,  etc.,  clusters  with 


T.  =  |  z  <l"...j"|t|l 
J  J  n>n" 


J 

.j’>  n 

i=l 


X 


i* 


Operator  T 


(II--) 


The  matrix  elements  of  T  are  the  unknown  numbers  which  define  the  formal 
operator  t.  The  number  j  corresponds  to  the  number  of  excitations: 

T  -  T,  +  T2  +  T3  +  T4  +  ...  ( 1 1 1 - 3 ) 

T-| :  1-electron  excitation  operator  (only  single  excitations) 

T2:  2-electron  excitation  operator  (only  double  excitations) 

T3:  3-electron  excitation  operator  (only  triple  excitations) 

T^:  4-electron  excitation  operator  (only  quadruple  excitations) 

To  make  some  connection  between  the  coupled  cluster  wave  function  of 
equation  (II-l)  and  more  conventional  Cl  (configuration  interaction)  or 
MBPT  (many-body  perturbation  theory)  expressions  for  |  ¥>,  expand  the 
right  side  of  equation  (II-l)  and  collect  terms  of  common  excitation  level: 

eT  I*  >  -  [1  +  (T,  +  T2  +  T3  +  T4  +...)  +  |,(T,  +  T2  +  Tj  +  T4  t...)2 

+  II(Tl  +  T2  +  T3  +  T4  tT!<T,  +  T2  +  T3  +  T4  +--->4] 


V  > 


( 1 1-4) 


By  grouping  the  terms  of  a  given  excitation  level  together,  we  see  that 
the  coupled  cluster  wave  function  can  be  rewritten  as 

eT  |  v  >  =  (1  +  C]  +  C2  +  C3  +  C4)|  v  >  ( 1 1-5) 

where  the  Cl  operators  C^ ,  C2>  C3>  C^,  ...  are 


•  *  o  . 


L4  =  4  Tl'l  +  1\  2 


1  3  21  ]2  ,etc. 


(H-6) 


The  method  has  been  coded  up  and  tested  for 
eT  =  t  +  T,  *  T2  +  +  T4  +  ^T,2  +  ll22  +  T,T2  + 

13  12  14 

+  3!T1  +  2T1  T2  +  4!T1  ^II_7) 

|  Sample  test  runs  have  been  made  on  H90  and  H.  using  all  of  the  above 

!  terms .  c  4 

2.  Coupled  cluster  method  (based  on  multi  reference  determinants) 

Since  both  our  preliminary  and  large  scale  ab-initio  configuration 
interaction  calculations  on  CH3NO2  indicated  a  2-configuration  ground 
state,  even  at  equilibrium  geometry,  and  even  more  configurations  along 
the  C-NO2  dissociation  pathway,  it  has  been  apparent  to  us  that  multi- 
reference  determinants  would  be  necessary  for  coupled  cluster  or  many- 
body  perturbation  theory  calculations  on  nitro  compounds. 

We  thus  initiated  as  a  first  step  the  derivation  of  necessary  formal 
equations  for  coupled  cluster  theory  based  on  multi  reference  determinants 
using 


T  =  T1  +  T2  (II-8) 

eT  =  1  +  (T-,  +  T2)  +  +  |t22  +  T-jT2  ( 1 1 -9) 

The  basic  premise  was  that  the  multi  reference  determinants  themselves 
already  contained  a  sizable  contribution  from  the  important  excitations. 
This  is  especially  true  for  the  way  we  obtain  our  multi  reference  deter¬ 
minant  which  is  from  ab-initio  configuration  interaction  calculations 
which  are  themselves  based  on  multi  reference  single  and  double  excitations 
from  the  important  multi  reference  determinants  from  a  smaller  preliminary 
Cl  calculation. 

The  equations  were  derived  at  the  Johns  Hopkins  University  by 
Professor  Laforgue  and  the  method  as  defined  by  equation  (1 1-9)  has 
been  coded  by  Dr.  Roszak,  debugged  and  is  now  undergoing  testing. 
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3.  Coupled  cluster  method  (repetitive) 

For  a  repetitive  couple  cluster  method,  e^  (using  only  lower  T . ' s ) 
would  be  applied  repetively  to  get  terms  to  any  desired  degree  of  J 
excitation.  The  advantage  of  such  a  repetitive  method  is  that  for  higher 
orders  of  excitation  one  does  not  expand  the  number  of  equations  that 
have  to  be  soved  at  the  same  time. 

As  the  first  step  in  implementing  this  procedure,  the  equations  were 
derived  (by  Professor  Laforgue  at  the  Johns  Hopkins  University)  using 

T  =  Tj  +  T2  (11-10 

eT  =  1  +  (T1  +  T2)  +  +  ^T22  +  T-|  T2  (11-11 

as  above. 


This  repetitive  method  as  defined  by  equation  (11-11)  has  been 
coded  by  Dr.  Roszak  for  a  single  reference  configuration,  debugged  and  is 
now  undergoing  testing.  The  repetitive  method  is  now  also  being 
extended  to  multi  reference  determinants. 
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C.  RDX  preliminary  results 

The  coupled  cluster  method  has  been  applied  to  RDX 

Hs/H 

02N— N  N— N02 

HsL  Lh 
j/  \n/  XH 
Ii02 


H- 


•# 


with  two  different  choices  of  starting  orbitals:  localized  and  con¬ 
ventional  canonical  delocalized  molecular  orbitals. 


1 .  Localized  orbitals 

For  preliminary  testing  purposes,  the  canonical  orbitals  from  an 
ab-initio  MODPOT/VRDDO  calculation  (66  basis  functions,  66  molecular 
orbitals  of  which  44  were  occupied  closed  shell  molecular  orbitals) 
were  localized  following  the  criterion  of  Bays.  From  the  66  space 
orbitals  (132  spin  orbitals),  a  subset  of  14  spin  orbitals,  of  which 
eight  were  occupied,  was  selected.  These  occupied  orbitals  were  local¬ 
ized  mainly  in  the  N-N  region.  From  this  set  of  orbitals  24  unique 
monoexcited  (T^),  216  di excited  configurations  (T2),  912  tri excited 
(T3)  and  1872  (t4)  configurations  can  be  generated. 

Even  using  just  the  linear  approximations,  the  number  of  simultan¬ 
eous  equations  to  be  solved  using  only 

eTi+T2  (11-12) 

Nx  =  24;  N2  =  216;  Nx  +  N2  =  240 

(  Nj  +  N2  )(  Nj  +  N2  )  =  57600  =  (240)2  (11-13) 

is  57,600,  which  is  the  total  number  of  elements  in  the  matrix. 

The  number  of  simultaneous  equations  to  be  solved  using 
gT i+T2+T  3+T4 

Nj  -  24;  N2  =  216;  N3  =  912;  N4  =  1872 

(Nx  +  N2  +  N3  +  N4)  =  3024 

(Hi  +  N2  +  N3  +  N4)  =9,144,576 
is  9,144,576. 

Within  the  computer  capacity  of  the  computer  on  which  these  calcu¬ 
lations  are  being  run  currently  (fixed  core  with  no  virtual  paging  and 
an  operating  limit  of  no  more  than  131072  (=  217)  words  in  any  one  ar¬ 
ray),  It  was  not  logistically  feasible  to  solve  9,144,576  simultaneous 
non-linear  equations.  We  analyzed  the  problem  and  a  program  written 
to  solve  a  system  of  equations  of  that  size,  which  would  have  to  go  in 
and  out  of  computer  core  memory  and  would  be  a  factor  of  the  order  of 
50  times  slower  than  an  in-core  program. 


(11-14) 

(11-15) 


Thus,  for  the  purposes  of  this  preliminary  testing,  we  used  the 
linear  approximations. 

e7  =  1  +  Ti  +  T2. 

Solutions  to  the  following  two  formulations  were  obtained: 

a.  Single  reference  function 
ECorr  =  -°-0105592  a-u- 

b.  Repetitive  correlation  approach 
E  =  -0.010600  a.u. 

[It  is  knowfrthat  the  linear  coupled  cluster  approach  may  overestimate 
the  correlation  energy  and  use  of  this  particular  linear  coupled  cluster 
approach  for  the  single  reference  was  just  for  comparison  purposes  to 
the  repeated  correlation  approach.] 

The  repetitive  method  is  being  derived  and  coded  at  present  for 
multi  reference  determinant  functions. 

2.  Canonical  orbitals 

A  large  scale  ab-initio  Cl  calculation  (using  multi  reference  deter¬ 
minants)  had  been  run  on  an  RDX  using  the  same  ab-initio  M0DP0T/VRDD0 
SCF  calculation  as  in  part  a.  Fourteen  space  orbitals  had  been  used 
for  that  Cl  calculation  (eight  of  which  were  occupied,  and  all  of  which 
had  appreciable  density  in  the  N-N  region)  allowing  all  single  and  double 
excitations  from  the  multi  reference  determinants  picked  from  a  prelimin¬ 
ary  smaller  Cl  calculation. 

The  two  dominant  configurations  results  from  that  multi  reference 
Cl  calculation  were  used  as  the  multi  reference  determinant  function  for 
this  coupled  cluster  calculation. 

Again  the  choice  made  for  these  coupled  cluster  calculations  was  the 
use  of  fourteen  (14)  molecular  spin  orbitals,  eight  (8)  of  which  were 
occupied.  These  spin  orbitals  had  appreciable  density  in  the  N-N  region. 
Again,  for  the  purposes  of  this  preliminary  testing  we  used  the  linear 
approximation  Z 

e  =  1  +  T!  +  T2 

a.  Single  reference  function: 

=  -0.0285416  a.u. 

corr 

(Since  the  canonical  orbitals  used  here  were  not  localized  as 
in  part  a.  but  only  had  appreciable  density  in  the  N-N  region  the  coupled 
cluster  correlation  energy  is  somewhat  larger  than  in  part  a.  where  the 
same  coupled  cluster  correlation  operator  and  a  single  reference  function 
were  used). 

b.  Mul ticonfiguration  reference  (ground  state  +  one  diexcited 
configuration  -  the  only  significant  terms  in  the  Cl  wave 
function) . 

Ecorr  =  -0-02854153  a.u. 
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Note:  In  this  case,  part  of  the  E  is  absorbed  into  the 
multi  configuration  SCF  energy.  co 


We  are  planning  to  rerun  the  ab-initio  Cl  calculations  using 
localized  orbitals  by  the  Boys  criterion.  We  will  then  carry  out  coupled 
cluster  calculations  on  the  Cl  multi  reference  determinant  wave  function. 

This  will  enable  a  comparison  of  single  reference  and  multi  reference 
coupled  cluster  based  on  localized  orbitals  and  also  a  comparable  comparison 
between  single  reference  and  repetitive  coupled  cluster  based  on  either 
localized  or  delocalized  molecular  orbitals. 

D.  Further  Investigations 

We  will  be  carrying  out  extensive  testing  of  the  repetitive  approach 
and  the  multi  reference  determinant  approach  to  coupled  cluster  theory, 
and  comparison  of  these  to  the  results  arising  from  testing  including  the 
various  terms  in  the  conventional  coupled  cluster  approach.  The 
objective  is  to  develop  a  tractable  desirable  computational  alternative 
to  the  brute  force  generation  and  solving  of  the  rapidly  ballooning  of 
millions  upon  millions  of  simultaneous  equations  to  be  solved  as  the 
systems  of  interest  become  larger  and  larger. 


w 

I 
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III.  Ab-initio  quantum  chemical  calculations  and  electrostatic  molecular 
potential  contour  maps  for  cationic  polymerization 

We  have  continued  optimizing  geometries  and  generating  electro¬ 
static  molecular  potential  contour  (EMPC)  maps  for  new  monomers  of 
interest  experimentally.  Our  previous  results  had  illustrated  vividly 
that  these  EMPC  maps  were  capable  of  predicting  reliably  the  basicity 
of  even  hypothetical  monomers  prior  to  their  synthesis.  The  initiation 
step  in  the  cationic  polymerization  of  cyclic  ethers  is  governed  by  the 
basicity  of  these  monomers  and  our  EMPC  maps  predicted  reliably  the 
rank  order  of  polymerization  and  even  optimal  copolymer  candidates. 

A.  Three-dimensional  electrostatic  molecular  potential  contour  (EMPC) 
maps  and  image  matching 

We  have  also  been  working  to  develop  methodology  for  more  quanti¬ 
tative  predictions  of  relative  polymerization  ratios.  This  is  complicated 
by  the  fact  that  there  are  two  steps  in  the  polymerization:  initiation-- 
proportional  to  basicity  and  propagation--considered  by  experimental 
polymer  chemists  to  be  influenced  by  the  electrophilicity  of  the  propa¬ 
gating  ring  opened  carbocation.  We  had  already  shown  that  these  two  steps 
are  influenced  in  opposite  and  competing  directions  by  the  character  of 
the  substituent  groups.  We  had  also  shown  that  initiation  was  a  dominant 
step.  Thus  we  have  developed  and  implemented  methods  to  calculate  the 
volume,  shape  and  accessibility  of  isopotential  EMPC  contours.  We  are  now 
carrying  out  such  calculations  to  compare  against  experimental  data.  We 
are  also  developing  a  strategy  for  comparing  el ectrophili cities  for  the 
propagation  step. 

A  further  objective  in  this  area  of  the  research  is  to  develop 
and  test  methods  to  predict  quantitatively  relative  polymerization  ratios. 
This  is  complicated  because  experimentally  measured  relative  polymerization 
ratios  are  governed  both  by  the  initiation  step  (influenced  by  the  basicity 
of  the  cyclic  ether  as  exemplified  by  its  EMPC  maps  and  supplementarily 
by  the  steric  accessibility  for  cationic  attack  at  the  oxygen)  and  by  the 
electrophilicity  of  the  propagating  carbocation. 

As  the  first  step  we  are  calculating  the  volume,  shape  and  accessibility 
of  some  of  the  prototype  exotic  energetically  substituted  oxetanes  such 
as  oxetane,  3-fluoro-3-nitrooxetane  and  3,3-dinitrooxetane.  The  relative 
volumes  and  other  quantities  are  being  calculated  within  various  specified 
EMPC  isopotential  contours  since  this  is  a  totally  new  theoretical  approach 
and  there  may  be  depending  on  particular  values  of  EMPC  contours. 


1.  Methodology 


An  electrostatic  molecular  potential  contour  map  is  the  set  of  points 


where  p(r' )  is  the  electron  density  at  r'  (obtained  from  our  ab-initio  MODPOT/ 
VRDDO/MERGE  SCF  programs),  is  the  charge  on  nucleus  a.  For  a  typical 

molecule,  V(r)  is  calculated  for  about  30,000  different  values  of  r  distri¬ 
buted  over  a  uniform  three-dimensjonal  grid.  V(r)  is  the  potential  a  posi¬ 
tive  test  charge  experiences  at  r.  Thus,  negative  values  of  V(r)  indicate 
electrophilic  regions,  and  positive  values  show  nucleophilic  regions.  We 
will  be  calculating  these  quantities  for  a  number  of  the  other  substituted 
oxetanes . 

We  are  limited  in  the  experimental  data  for  relative  polymerization 
ratios  against  which  to  compare.  We  will  be  carrying  out  such  calculations 
for  the  oxetanes  for  which  the  experiments  have  been  run. 

There  are  several  different  approaches  for  computer  image  matching  of 
EMPC  maps  for  quantitative  (or  semi-quantitative)  prediction  of  relative 
polymerization  ratios. 

a.  Volume 

Match  the  volumes  for  various  isopotential  EMPC  contours  of  differently 
substituted  oxetanes. 

b.  Shape 

Compare  the  radii  at  which  various  isopotential  EMPC  contours  of  dif¬ 
ferently  substituted  oxetanes  intersect  a  common  set  of  vectors  directed 
outward  from  designated  points  on  the  molecule  (for  example,  fron  the  ring 
oxygen  atom). 


c.  Accessibility 

The  regions  of  clear  access  to  the  various  EMPC  contours  around  the 
ring  oxygen  atom  without  interference  from  other  parts  of  the  molecule. 


EMPC  points  are  generated  in  a  rectangular  grid  spanning  each  plane 
The  location  and  the  number  of  the  planes  should  be  such  as  to  give 
a  complete,  closed  picture  of  the  contour  of  interest  (else  the 
contour  will  be  “clipped"). 


(1).  Determine  volume  of  a 
typical  contour,  say 
-20  Kcal/mole. 


(2).  Draw  the  vectors  indicated,  at  least  two  in  each  coordinate  plane. 
These  are  tangents  from  a  selected  point  (01,  in  the  example)  to  the 
contours. 


1 


(3).  This  set  is  a  descriptor  of  the  contour  under  consideration.  It  is 
possible  to  augment  this  descriptor  set  by  additional  vectors 


2.  Preliminary  Results 


a.  Volume 

The  following  results  indicate  that  while  the  volume  of  the  -20  kcal 
contour  OXET  >  FNOX  >  DNOX,  the  averaged  coordinates  of  the  centroids 
of  volumes  are  not  very  different. 

VOLUME  OF -20  KCAL/MOLE  CONTOUR 
OXET  21.90  A3 

FNOX  3.54  A3 

DNOX  1.26  A3 

AVERAGED  COORDINATES  OF  THE  VOLUME 
OXET  0,  0,  -2.37 

FNOX  0,  0.01,  -2.17 

DNOX  0,  0,  -2.17 

Thus  another  approach  was  tried;  namely,  a  measure  of  the  shape  and 
accessibility  of  EMPC  isopotential  contours. 
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B.  Molecules  Calculated 

We  are  still  carrying  out  ab-initio  MODPOT/VRDDO/MERGE  calculations 
for  geometry  optimizations  and  generating  EMPC  maps  for  others  of  the 
large  number  of  new  substituted  oxetanes  of  interest  for  cationic  poly¬ 
merization  which  Jerry  Manzer  had  discussed  with  us  on  his  visit  to  us  at 
the  Johns  Hopkins  University. 

What  makes  the  quantum  chemical  calculations  so  involved  for  most  of 
the  molecules  Manzer  mentioned  is  the  number  of  multiple  maxima  and  minima 
for  their  energies  as  a  function  of  geometry. 

Since  so  many  geometries  have  to  be  calculated  for  these  molecules, 
minor  differences  in  integral  prescreening  and  cut-off  values  became 
significant.  We  carried  out  extensive  calculations  with  different  integral 
threshholds  and  have  established  the  appropriate  criteria. 

Many  of  the  molecules  are  difficultly  convergent  in  the  SCF  procedures. 
This  difficult  convergence  had  also  affected  our  previous  calculations  on 
the  calculations  for  proton  attack  on  the  cyclic  ethers.  We  implemented 
new  SCF  routines  with  special  convergence  enhancing  procedures  which  have 
overcome  this  problem. 

We  have  carried  out  further  calculations  with  these  above  improvements 
on  a  number  of  the  cyclic  ethers  we  previously  investigated  to  assess  the 
importance  of  such  techniques.  These  results  will  be  incorporated  into  the 
technical  manuscripts  which  are  being  written  or  prepared. 

Several  manuscripts  have  been  published.  In  our  presentations  at 
meetings  we  display  the  three-dimensional  EMPC  maps  in  a  combination  of 
various  colors  which  give  a  vivid  portrayal.  It  is  difficult  to  portray 
three-dimensional  EMPC  maps  in  black  and  white  for  publication.  The 
only  feasible  way  to  portray  these  in  black  and  white  is  to  display  each 
three-dimensional  EMPC  isopotential  contour  seperately,  which  is  less 
than  desirable.  We  have  been  talking  to  publishers  of  several  leading 
journals  who  have  expressed  possible  interest  in  looking  into  publishing 
our  EMPC  maps  in  color  provided  we  can  provide  appropriate  color  separates 
and  providing  their  costs  for  such  color  pages  is  not  too  excessive. 

During  the  past  year,  presentations  of  this  research  have  been  given 
at  a  number  of  national  and  international  meetings  and  at  other  univer¬ 
sities  and  research  institutions. 

At  the  invitation  of  Dr.  Eli  Pearce,  Chairman  of  the  1983  Gordon 
Conference  on  Polymers,  Dr.  Kaufman  presented  the  results  of  our  research 
on  quantum  chemical  investigations  of  the  mechanism  of  cationic  polymer¬ 
ization  at  that  Gordon  Conference,  "Quantum  Chemical  Calculations  on  the 
Mechanism  of  Cationic  Polymerization,"  July  1983.  Our  theoretical  re¬ 
search  is  considered  to  be  both  pioneering  and  of  forefront  significance 
in  the  area  of  polymerization.  There  was  a  great  deal  of  interest  expressed 
from  a  number  of  industrial  companies  who  manufacture  polymers. 


By  invitation.  Dr.  Kaufman  also  presented  an  invited  plenary  lecture 
on  the  results  on  polymers,  "Ab-Initio  MODPOT/VRDDO/MERGE  Calculations 
(Intra-  and  Intermolecular)  for  Polymers  and  Solids,"  at  the  NATO  Advanced 
Study  Institute  on  Quantum  Chemistry  of  Polymers:  Solid  State  Aspects, 
Braunlage/Herz,  West  Germany,  July  1983. 

By  invitation.  Dr.  Kaufman  also  presented  these  results: 

"Three-Dimensional  Electrostatic  Molecular  Potential  Contour  Maps 

A.  Stereoelectronic  Requisites  for  Biomedical  Molecules 

B.  Cationic  Polymerization,"  8th  Canadian  International 

Theoretical  Chemistry  Symposium,  Halifax,  Canada,  August  1983 

Dr.  Kaufman  will  also  be  presenting  these  polymer  results  at  the  BF 
Goodrich  Company  in  October  1983. 


Credit  for  the  calculations  on  the  polymers  and  the  computer  image 
matching  is  due  to  Dr.  P.  C.  Hariharan. 
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IV.  Prediction  of  Crystal  Density 
A.  Background 

Detonation  pressures  and  detonation  velocities  are  influenced  by  crystal 
densities. 


To  be  able  to  predict  crystal  densities  even  for  completely  hypothetical 
compounds  is  of  interest  to  ONR.  While  there  are  empirical  methods  to  pre¬ 
dict  crystal  densities,  these  methods  depend  only  on  the  number  and  type 
of  fragments  or  groups.  These  empirical  methods  are  incapable  of  predicting 
differences  in  crystal  densities  between  various  position  isomers  of  the 
same  molecule.  These  empirical  methods  for  crystal  densities  are  also 
utterly  incapable  of  predicting  crystal  structure  arrangements. 

Thus,  we  embarked  several  years  ago  on  a  project  to  calculate  optimal 
crystal  packing  and  crystal  structure  parameters  based  on  potential 
functions  from  energy  partitioned  ab-initio  intermolecular  SCF  calculations 
plus  calculations  of  dispersion  energy  contributions. 

Our  approach  aimed  at  evaluating  the  intermolecular  interactions  is 
based  on  nonempirical  ab-initio  calculations  for  smaller  molecular  aggregates 
(monomers,  dimers,  trimers,  etc.),  partitioning  the  total  SCF  interaction 
energies  into  the  different  components  and  then  fitting  these  components 
individually  to  functional  forms  or,  when  necessary,  recalculating  or  estim¬ 
ating  explicitly  for  certain  interaction  components  for  each  different  unit 
cell  dimension  change  plus  calculated  estimates  of  dispersion  energy. 

The  general  method  was  described  in  detail  in  our  1982  Annual  Report  on 
this  project. 

The  main  components  of  the  two-body  interaction  energy  are 


aE  *  +  E^  +  E^ 

A  tEL  lEX  tIND  tDISP 


(IV-1) 


The  first  three  contributions  to  equation  IV-1  can  be  obtained  from  the 
decomposition  of  the  variational  E<.rn  interaction  energy  corrected  for  basis 
set  superposition  error  (BSSE). 

&ESCF  =  EAB  "  EA(B)  "  EB(A) 

*  c(l)  +  r(l)  +  c(R)  +  f(R) 

=  LEL  lEX  iND.LE  tIND,CT 

where 

EAB*  ea(B) ’  EB( A) *  denote  the  SCF  total  energies  of  interacting 
dimer,  isolated  monomers  A  and  B  (evaluated  in  the  dimer  basis  set)  respectively 


and  where 


Eel  -  the  electrostatic  contribution  is  further  partitioned 

.  E0>  +  Etl) 

EL.MTP  ‘El, PEN 


’.yoo- ..  .. 


tiY'  -  the  exchange  contribution 


E 


(2) 

IND 


-  the  induction  contribution  is  further  partitioned 


=  E{2) 
tIND,LE 


+  E 


(2) 

IND.CT 


e(2) 

DISP  -  the  dispersion  energy  contribution. 


In  addition,  the  induction  component  contains  also  remaining  (R)  higher 
order  correction  (R)  and  it  is  decomposed  into  an  additive  long  range 
classical  component  E[^  ^  arising  from  interaction  of  local  excitations 
(LE)  with  static  charge  1  L  distributions  of  other  molecules  and  a  short 
range  EV^  qj  corresponding  to  interaction  of  ionic  charge-transfer  config¬ 
urations  ’  (CT)  with  static  charge  distributions. 


B.  Developments  this  Past  Year 
1.  Programs 

a.  MOLASYS- INTER 

(1) .  Calculation  of  atomic  quadruples  for  use  in  calculating 

E^MjP  and  in  calculating  the  electric  field  has  been 
implemented,  tested  and  used; 

(2) .  An  efficient  initial  guess  has  been  added  (which  amelio- 

ates  previous  problems  with  non-convergence  for  some 
systems) ; 

(3) .  Classical  perturbation  and  variation  perturbation  ex¬ 

pressions  for  calculating  dispersion  energy  have  been 
implemented  and  tested. 

b.  CRYSTAL 

(1) .  Use  of  atomic  quadrupoles  in  calculating  E^MTP  has 

been  implemented,  tested  and  used; 

(2) .  Calculation  of  the  electric  field  up  through  quadrupole 

terms  at  the  molecules  in  the  reference  cell  from  all 
molecules  in  all  surrounding  unit  cells  has  been  imple¬ 
mented,  tested  and  used; 

(3) .  Calculation  of  two-body  and  three-body  induction  terms 

from  the  calculated  electric  field  and  atomic  polariza¬ 
bilities  has  been  implemented,  tested  and  used; 

(4) .  Stepwise  increase  of  crystal  cluster  size  has  been 

implemented  and  tested  in  order  to  investigate  the  con¬ 
vergence  of  the  results  with  the  increase  of  crystal 
cluster  size. 
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2.  Ab-Initio  atom  class-atom  class  potentials 

a.  Nature  of  anisotropy  of  non-empirical  ab-initio  atom-atom 
potentials 

Most  of  the  currently  employed  intermolecular  atom-atom  potential  func¬ 
tions  (either  empirical  or  theoretical)  are  assumed  to  be  isotropic,  and 
very  little  is  known  about  the  nature  of  their  anistropy  on  an  atomic  level. 
This  situation  stronly  limits  the  applicability  of  available  potentials  for 
studying  properties  of  condensed  matter,  etc. 

Based  on  our  ab-initio  interaction  energy  components  calculated  for 
the  model  CO  dimer,  we  analyzed  the  nature  of  the  anistropic  contribution 
expressed  in  the  form  of  atom-atom  potentials.  For  this  purpose,  all  the 
main  contributions  of  the  intermolecular  interaction  energy  have  been  re¬ 
presented  in  the  form  of  atom-atom  Q(a)  R“J2  potentials  for  different 
angular  orientations  (a)  of  the  CO  dimer.  (While  the  R^J2  form  is  not 
optimal  for  all  components,  this  gives  us  a  chance  to  compare  directly  the 
anisotropy  of  all  contributions  simultaneously.)  Results  (Figure  1.)  indi¬ 
cate  that  the  observed  anisotropy,  i.e.,  a  dependences  of  Q,  is  almost 
entirely  due  to  the  electrostatic  multi  pole  terms,  whereas  the  exchange, 
electrostatic  penetration,  induction  and  dispersion  contributions  could  be 
successfully  represented  by  isotropic  atom-atom  potentials. 


This  reinforces  our  approach  where  we  calculate  the  strongly  anisotropic 
electric  multipole  terms,  Ej.')  directly  from  the  cumulative  atomic  mul- 

,  EL)  mr 

tipole  moments  derived  from  the  molecular  wave  function. 

The  other  terms  Eg]^»  and  CT,  we  fit  by  the  following 

isotropic  potential 

E  =  (b  +  yR-1 )exp(-$R) 

Besides  the  electrostatic  multipole  term,  the  next  most  anisotropic 
contribution  seems  to  be  dispersion  interaction,  but  to  a  much  smaller  degree. 
This  term  is  better  fitted  by  Qg(a)  R~^  potential  where  Qg  also  displays 

some  a  dependence. 


b.  Non-empirical  atom-atom  potentials  for  main  components  of 
intermolecular  interaction  energies 


The  aim  of  this  research  was  to  obtain  non-empirical  atom-atom  potential 
for  each  main  term  of  intermolecular  interaction  energies  by  separately 
fitting  ab-initio  SCF  interaction  energy  partitioned  energy  components 
evaluated  for  a  number  of  model  dimers  (symmetric  or  mixed)  within  our 
energy  decomposition  scheme  corrected  for  basis  set  superposition  error 
(BSSE).  The  calculations  have  been  performed  with  well-balanced  minimal 
valence  3S3P  MODPOT  basis  sets  with  ab-initio  effective  core  model  potentials 
known  to  yield  interaction  energies  very  close  to  results  obtained  with 
extended  basis  sets  at  much  less  computer  time.  The  non-transferable 
^IND  CT^  ant*  stron9ly  anistropic  (E^  MTP^  liave  *3een  calculat;ed  explicitly 

each  time  from  the  cumulative  atomic  multi  pole  moments  (CAMM),  up  through 
quadrupole  terms,  obtained  directly  from  the  monomer  wave  functions.  The 

short  range  exchange  (E^),  ctiar9e  transfer  (e|^  ct)  and  electrostatic 

penetration  (E^PEN)»  have  been  represented  by  simplified  potentials  in 

the  form  of  (g  +  yR-l)exp  ( -6R) ,  fitting  corresponding  results  for  dimers 
(symmetric  or  mixed)  formed  by  HF,  HoO,  NH3,  CH4,  CO,  C02.  These 
particular  molecules  were  chosen  as  the  most  appropriate  from  which  to 
obtain  potentials  for  the  first  row  atoms  (C,N,0,F)  and  H. 

The  polarizability  type  terms  oa  and  used  in  calculation  of  the 


EIND  LE  from  the  calcu1ated  value  of  the  electric  field  and  a  were  obtained 

by  equating  E^  lc  directly  from  the  energy-partitioned  SCF  wave  function 
inu,Lt  (2)AB 

to  the  expression  for  EJ^  ^  calculated  using  the  value  of  the  electric 
field  calculated  within  the  atomic  multi  pole  expansions. 


.(2)  AB 
■IND.LE 
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c.  Transferability  of  short  range  potentials  for  nitro  compounds 

From  a  large  number  of  intermolecular  ab-initio  SCF  calculations  on 
(CH3NO2)?  we  fitted  the  short  range  ab-initio  atom  class-atom  class 
potentials  for  00,  OH,  NO  and  HH  contacts.  Similarly,  from  a  large  number 
of  intermolecular  ab-initio  SCF  calculations  on  the  nitrobenzene,  CgHd^, 
dimer  we  fitted  ab-initio  atom  class-atom  class  potentials  for  the  U0, 
OH,  NO  and  HH  contacts. 


FIGURE  2. 


TRANSFERABILITY  OF  ATOM -ATOM  POTENTIALS 
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C6H5N02  potential  for  O-H 
. CH3N0?  potential  for  O—H 


As  can  be  seen  from  Figure  2,  there  is  very  close  agreement  in  the 

:0>  cO)  c(2) 


fitted  ab-initio  short  range  0 — H  potentials  E 


EX 


cv 1 /  and  C\ 

"EL.PEN  ana  tIND,CT  from 
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C.  Results 

Using  the  short  range  potentials  fitted  from  ab-initio  intermolecular 

calculations  on  the  (CHJK^g  dimer  and  explicit  calculation  of  Eg^Mjp  and 

the  electric  field  for  calculation  of  ic  from  the  ab-initio  SCF  cal- 

INDjLE 

culation  on  RDX 


H 

H 


>C 


Hv  H 


N  -  NOr 


.CC, 


we  calculated  the  unit  cell  parameters  of  RDX  to  -1.0%(a),  1.5%(b),  -0.6%(c) 
and  the  crystal  density  to  0.1%  deviation  from  the  experimentally  measured 
crystal  structure. 
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Major  credit  is  due  to  Dr.  W.  Andrzej  Sokalski  for  the  improvements 
to  the  CRYSTAL  program  and  the  potentials.  Credit  is  also  due  to  Dr.  Alfred 
H.  Lowrey,  Naval  Research  Laboratory,  who  is  working  with  us  collaboratively 
on  the  project  as  a  fellow-by-courtesy  of  the  Johns  Hopkins  University. 

Dr.  Szczepan  Roszak  and  Dr.  John  M.  Blaisdell  also  contributed  to  this 
project. 


V.  Ab-initio  calculations  on  polynitropolyhedranes  and  polynitropolyaza- 

polyhedranes 

We  have  continued  to  carry  out  ab-initio  calculations  on  polynitropoly¬ 
hedranes  and  polynitropolyazapolyhedranes .  The  calculated  quantum  chemical 
indices  such  as  total  overlap  populations  (TOP's)  between  atoms  (skeletal  and 
skeletal -NO2)  as  a  function  of  nitration  indicate  whether  and  where  nitration 
affects  inherent  bond  strengths.  Also,  from  these  wave  function  we  have 
been  generating  electrostatic  molecular  potential  contour  (EMPC)  maps  around 
these  polynitropolyhedranes  and  polynitropolyazapolyhedranes.  These  EMPC  maps 
appear  to  indicate  that  the  character  of  the  nitro  groups  may  change  as  a 
function  of  the  degree  of  nitration.  Also  these  EMPC  maps  may  prove  useful  in 
understanding  how  these  molecules  will  pack  in  crystals. 

Molecules  which  have  a  number  of  nitro  groups  around  their  peripheries 
could  scarcely  be  expected  to  approach  each  other  closer  than  a  zero  iso¬ 
potential  contour  if  the  signs  of  the  EMPC  contour  of  the  adjoining  mole¬ 
cules  are  the  same.  On  the  other  hand,  if  the  molecules  are  surrounded  by 
by  molecules  with  both  positive  and  negative  EMPC  contours,  the  molecules 
may  try  to  orient  themselves  in  a  crystal  having  regions  of  positive  EMPC 
contours  around  one  molecule  adjoining  regions  of  negative  EMPC  contours 
around  the  other  molecule. 

A.  Electrostatic  molecular  potential  contour  (EMPC)  maps  around 
nitrocubanes 

As  an  indication  that  the  character  of  the  nitro  groups  might  be  changing 
as  a  function  of  the  degree  of  nitration,  the  EMPC  maps  for  1 ,3,5,7-tetranitro- 
cubane  and  for  octanitrocubane  are  presented  on  the  next  pages  (Figures  4  and  5). 

The  small,  very  dark  2,000  kcal/mole  isopotential  contours  are  "black" 
and  denote  the  positions  of  the  nuclei.  The  next  darker  -10  kcal/mole  iso¬ 
potential  contours  are  located  around  the  nitro  groups.  However,  in  the 
1 ,3,5,7-tetranitrocubane,  these  -10  kcal/mole  EMPC  regions  are  rather 
delocalized  around  the  entire  N0«  group.  In  octanitrocubane  (which  is  in 
exactly  the  same  orientation  as  The  1 ,3,5,7-tetranitrocubane) ,  the  -10 
kcal/mole  contours  in  the  regions  of  the  N0p  groups  are  much  more  sharply 
localized  over  the  individual  atoms. 


B.  Polynitroadamantanes 

We  have  carried  out  ab-initio  M0DP0T/VRDD0  calculations  on  other 
polynitroadamantanes  whose  crystal  structure  had  been  determined 
experimentally,  and  also  from  these  calculated  wave  functions  generated 
the  three-dimensional  isopotential  EMPC  maps. 

Among  the  compounds  investigated  were  tetranitroadamantane  and 
gens-dinitroadamantane. 

This  will  give  us  some  feeling  for  the  magnitudes  of  various  quantum 
indices  for  compounds  that  do  exist. 

Major  credit  for  these  calculations  is  due  to  Dr.  P.C.  Hariharan. 
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VI.  Ab-initio  configuration  interaction  calculations  on  the 
decomposition  pathway  of  nitromethane 

Ab-initio  configuration  interaction  calculations  using  multirefer¬ 
ence  determinants  have  been  carried  out  for  much  of  the  CH^-NO^  decomp¬ 
osition  pathway  of  nitromethane.  For  each  point  in  the  order  of  800,000 
configurations  were  generated.  Using  a  perturbation  procedure  the 
energy  contribution  of  each  configuration  was  calculated.  All  configur¬ 
ations  contributing  to  an  energy  higher  than  a  chosen  threshold  were 
included  in  the  Cl  wave  function  to  be  solved.  Then  an  extrapolation 
was  performed  to  take  into  account  the  energy  contributions  of  the  con¬ 
figurations  not  included  explicitly.  This  was  then  followed  bya  David¬ 
son  type  correction  for  size  consistency. 

A.  Details  of  configuration  interaction  calculations 

Reference  configurations  =19  (All  single  and  double  excitations  from 

these  19  reference  configurations) 

Electronic  states  explored  =  4 

Total  SAF's  generated  =  845,020  (at  equilibrium  C-N  distance  2.8  a.u.) 
Interaction  (selection)  threshold  =  0.000050  Hartree 
r(C-N)  a.u.  |  2.3  |  2.8  I  3.3  1  3.8  1  4.3  1  4.8  ! 


B.  Results 


There  are  several  significant  aspects  of  these  Cl  results  on  CH3N02. 

1.  Multiccnngurational  character  even  of  ground  electronic  state  at 
equilibrium  geometry 

First,  even  at  the  equilibrium  geometry  for  the  electronic  ground 
state  of  CH^NO^  there  is  a  significant  multi  configuration  character. 

This  is  due  to  a  double  excitation  localized  out-of-plane  on  the  NCL 
group. 

As  the  r^._N  is  increased,  even  additional  configurations  become  im¬ 
portant  (corresponding  to  breaking  the  C-N  a  bond)  in  addition  to  the 
multiconfigurational  character  localized  on  the  N0„  group  which  always 
remains  important.  L 

For  the  first  excited  and  higher  excited  states,  there  is  even  more 
mixing  and  even  stronger  multi  configuration  character  at  each  rr  N  dis¬ 
tance.  L-IN 

Thus,  even  for  the  ground  electronic  state  of  CH3N02  multi  reference 

determinants  will  be  necessary  for  Cl  calculations  or  for  many-body 
perturbation  theory  or  for  coupled  cluster  theory  calculations. 

2.  The  molecular  orbital  ordering  of  both  occupied  and  virtual 
orbitals  changes  between  3.3  and  3.8  a.u. 

We  denote  the  orbitals  described  as 

ac_N  bonding  in  the  C-N  bond  in  the  plane  of  the  N02  group 

a*c_N  anti  bonding  in  the  C-N  bond  in  the  plane  of  the  N02  group 


irn  n  occupied,  localized  primarily  on  0,0,  perpendicular  to  the 

ulu2  plane  of  the  N02  group  c 

u*Nn  n  unoccupied,  localized  primarily  on  N0,0,  perpendicular  to 
<ulu2  the  plane  of  the  N02  group  6 

HOMO  LUMP  2nd  LUMP 

^  02  1T*N0102  °*c-n 

°C-N  °*C-N  1T*N0102 


2nd  HOMO 

-  3.3  a.u.  a, 


=  3.8  a.u. 


'C-N 

ff0l02 


Major  credit  for  these  Cl  calculations  is  due  to  Dr.  Cary  Chabalowski. 


VII.  Ab-initio  configuration  interaction  calculations  on  RDX 

Preliminary  Cl  calculations  were  carried  out  on  RDX  at  equilibrium 
considering  the  molecular  orbitals  centered  primarily  on  a 

c^n-no2 

portion  of  the  skeleton  as  a  prototype  for  the  N-NO^  decomposition. 

Even  at  equilibrium  a  multi  reference  determinant  character  on  the  N0? 
group  was  found— similar  to  that  in  CH^C^.  Thus,  it  appears  to  be  L 

a  general  characteristic  of  all  aliphatic  nitroexplosives. 

A.  Details  of  the  calculation 


For  this  preliminary  test  Cl  calculation  on  RDX 


\  /H 

0,N— N  N— NO,, 


H 
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C  C 

W 

no2 


we  used  our  customary  ab-initio  MODPOT  basis  set  and  the  VRDDO  prescreen¬ 
ing  option  which  cuts  down  considerably  on  the  number  of  integrals  to  be 
calculated.  There  were  42  occupied  molecular  orbitals.  Of  these,  17 
occupied  molecular  orbitals  were  retained  as  well  as  15  virtual  orbitals 
for  the  configuration  interaction  calculation.  These  molecular  orbitals 
were  ones  with  their  dominant  contributions  in  a  particular  one  of  the 

C\ 

>N— NO, 

Z/  6 

regions  as  a  model  for  a  r 

> . NO. 

c  c 

decomposition  pathway.  The  ab-initio  Cl  calculation  was  a  multi  reference 
determinant  calculation  allowing  all  single  and  double  excitations  from 
19  reference  configurations.  The  total  number  of  SAF's  generated  was 
368,845.  The  energy  contribution  from  each  of  these  SAF's  was  estimated 
by  a  perturbation  procedure  and  all  SAF's  contributing  an  energy  more 
than  a  certain  threshold  were  retained  (3360  configurations)  for  the 
Cl  calculation.  Four  roots  were  obtained  corresponding  to  the  ground 
and  first  three  excited  singlet  configurations.  Then  the  energy  contri¬ 
butions  from  the  remaining  configurations  were  added  by  extrapolating. 
Following  that  a  Davidson  type  correction  was  added  for  size-consistency. 
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Even  in  the  ground  electronic  state  at  the  equilibrium  geometry  of 
RDX  there  were  contributions  from  the  doubly  excited  closed  and  open 
shell  configurations  localized  on  the  NO^  group  as  well  as  smaller  con¬ 
tributions  from  other  configurations  in  that 

C^N-N0o 
C"*  1 

region. 

A  Cl  calculation  on  the  RDX  molecule  taking  into  account  simultan¬ 
eously  the  three 

C^N-N0o 

2 

regions  would  be  expected  to  have  similar  contributions  from  all  three 
groups  simultaneously. 

We  will  be  investigating  further  Cl  calculations  on  the  decomposi¬ 
tion  of  RDX. 

Major  credit  for  these  Cl  calculations  is  due  to  Dr.  Cary  Chabalowski 
with  the  collaboration  of  Dr.  P.C.  Hariharan. 


«.  *\.  ‘. 


1  w.  T 1  c  , 
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VIII.  Conversion  of  our  programs  to  vector  supercomputers,  CYBER  205 

and  CRAY 

The  quantum  chemical  problems  we  are  investigating  involve  large 
scale  calculations  on  large  systems,  both  with  many  atoms  and/or  with 
many  configurations  or  simultaneous  equations  to  be  solved.  We  have 
gotten  access  to  CYBER  205  time  and  to  CRAY  time.  The  CYBER  205  has 
two  million  64  bit  words  in  real  memory  plus  virtual  memory.  The 
CRAY  has  one  million  64  bit  words  of  real  memory.  We  are  in  the  process  of 
converting  our  programs  to  both  of  these  machines. 

We  are  enthusiastic  about  computing  on  these  supercomputers  because 
they  will  enable  us  to  do  many  problems  we  are  investigating  in  the  ONR 
(and  ARO)  research  which  have  been  virtually  intractable  in  more  con¬ 
ventional  computers. 

However,  these  conversions  of  our  programs  to  the  CYBER  205  and 
to  the  CRAY  are  non-trivial  even  just  to  convert  the  programs  to  get 
them  running  even  before  vectorizing  the  programs  for  optimal  speed 
since: 

A.  Our  programs  are  large  and  to  run  them  in  conventional  machines 
we  had  to  make  extensive  use  of  moving  data  in  and  out  of  mass 
storage  which  is  now  only  partly  necessary  and  is  done  in  a 
different  manner  on  the  CRAY  and  in  a  yet  different  manner  on 
the  CYBER  205. 

B.  Our  programs  made  extensive  use  of  masking  and  bit  shift  oper¬ 
ations  which  are  done  in  a  different  manner  on  the  CRAY  and 

in  yet  a  different  manner  on  the  CYBER  205. 

C.  The  permissible  FORTRAN  on  the  computers  we  have  been  running 
(CDC  6600,  CDC  7600,  CYBER  175)  is  somewhat  different  than  the 
permissible  FORTRAN  on  the  CRAY  and  is  somewhat  different  than 
the  permissible  FORTRAN  on  the  CYBER  205.  (The  permissible 
FORTRAN  on  the  CRAY  and  on  the  CYBER  205  are  al so  di fferent  from 
one  another.) 

D.  The  CYBER  205  is  a  newer  computer  than  the  CRAY  and  hence  its 
operating  system  is  still  undergoing  developmental  changes 
which  sometimes  necessitates  changes  in  the  program. 

However,  we  have  already  sucessfully  converted  one  of  our  programs 
and  we  are  working  on  the  others.  We  will  then  work  on  proper  vector- 
ization  of  the  code  for  optimal  running  speed. 


Major  credit  is  due  to  Dr.  P.C.  Hariharan  in  this  endeavor. 


45 


IX.  Lectures  Presented  on  this  ONR  Research  (*denotes  invited  lecture) 

A.  Already  Presented,  October  1,  1962  -  September  30,  1963 
Dr.  Joyce  J.  Kaufman 

1.  At  National  or  International  Meetings 

*  "Ab-Initio  MODPOT/VRDDO/MERGE  Intermolecular  Calculations  for 
Crystals,"  Sanibel  International  Symposium  on  Molecular  and  Solid  State 
Theory  and  Computational  Quantum  Chemistry,  Palm  Coast,  Florida,  March  1983. 

*  "Ab-Initio  MODPOT/VRDDO/MERGE  SCF  and  Beyond  Hartree-Fock 
Calculations  on  Energetic  Compounds,"  Sanibel  International 
Symposium  on  Molecular  and  Solid  State  Theory  and  Computational  Quantum 

Chemistry,  Palm  Coast,  Florida,  Marcn  1983, 

"Ab-Initio  MODPOT/VRDDO/MERGE  Intermolecular  Calculations  to 
Derive  Atom  Class-Atom  Class  Potential  Functions,"  American 
Physical  Society  Spring  Meeting,  Baltimore,  Maryland,  April 
1983. 

"Prediction  of  Crystal  Densities  Using  Atom  Class-Atom  Class 
Potential  Functions  from  Energy  Partitioned  Ab-Initio  SCF 
Calculations,"  American  Physical  Society  Spring  Meeting,  Baltimore, 

Maryland,  April  1983. 

*  "Quantum  Chemical  Calculations  on  Mechanism  of  Cationic  Poly¬ 
merization,"  Gordon  Conference  on  Polymers,  New  London,  Conn., 

July  1983. 

*  "Ab-Initio  MODPOT/VRDDO/MERGE  Calculations  (Intra-  and  Inter¬ 
molecular)  for  Polymers  and  Solids,"  NATO  Advanced  Study  Institute 
on  Quantum  Chemistry  of  Polymers:  Solid  State  Aspects,  Braunlage/ 

Herz,  West  Germany,  July  1983. 

*  "Prediction  of  Crystal  Densities  Using  Ab-Initio  Potential 
Functions  from  Energy-Partitioned  Ab-Initio  M0DP0T/VRDD0  Inter¬ 
molecular  SCF  Calculations  Plus  Dispersion,"  8th  Canadian 
Theoretical  Chemistry  Symposium,  Halifax,  Canada  August  1983. 

"Three-Dimensional  Electrostatic  Molecular  Potential  Contour  Maps 

A.  Stereoelectronic  Requisites  for  Biomedical  Molecules 

B.  Cationic  Polymerization,"  8th  Canadian  International  Theoretical 
Chemistry  Symposium,  Halifax,  Canada,  August  1983. 
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IX.  Lectures  Presented  on  this  ONR  Research 


(*denotes  invited  lecture) 


A.  Already  Presented  (cont.) 


"Ab-Initio  MODPOT/VRDDO/MERGE  Intermolecular  Calculations  to 
Derive  Atom  Class-Atom  Class  Potential  Functions,"  American 
Chemical  Society  Fall  National  Meeting,  Washington,  D.C.,  Sep¬ 
tember  1983. 

"Prediction  of  Crystal  Densities  Using  Atom  Class-Atom  Class 
Potential  Functions  from  Energy  Partitioned  Ab-Initio  SCF  Calcu¬ 
lations,"  American  Chemical  Society  Fall  National  Meeting, 
Washington,  D.C.,  September  1983. 

"Ab-Initio  MODPOT/VRDDO/MERGE  SCF  Calculations  on  Energetic 
Compounds:  Nitrocubanes,"  American  Chemical  Society  Fall 
National  Meeting,  Washington,  D.C.,  September  1983. 

"Ab-Initio  Cl  Calculations  on  CH3NO?  Confirm  the  Mul ticonfigur- 
ational  Character  of  Both  Ground  ana  Electronically  Excited  States 
Even  at  Equilibrium  Geometries,"  American  Chemical  Society  Fall 
National  Meeting,  Washington,  D.C.,  September  1983. 


B.  Invited  Lectures  Presented  at  ONR  or  DOD  Meetings:  Joyce  J.  Kaufman 

"Ab-Initio  Quantum  Chemical  Calculations  on  Large  Polynitro-Polyhe- 
dranes;  Cl  Calculations  on  the  Dissociation  Pathway  of  CH3NO2; 

Prediction  of  Crystal  Densities,  ONR  Workshop  to  Review  Dense  Explosives 
Research  Program  and  Opportunities,  China  Lake,  California,  April  1983. 

"Quantitative  Measurement  from  Electrostatic  Molecular  Potential 
Contour  Maps  of  Oxetanes,"  Energetic  Polymer  Workshop,  Great  Oak, 
Maryland,  July  1983.  (Dr.  P.C.  Hariharan) 

"Ab-Initio  Calculations  on  Nitroexplosives  and  Prediction  of  Crystal 
Densities,"  Chemical  Decomposition  of  Energetic  Materials,  Great  Oak, 
Maryland,  August  1983. 
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X.  Papers  Published  on  This  ONR  Research,  October  1982  -  September  1983 
A.  Already  Published 

1.  W.A.  Sokalski,  S.  Roszak,  P.C.  Hariharan  and  Joyce  J.  Kaufman, 

"An  Improved  SCF  Interaction  Energy  Decomposition  Scheme  Corrected 
for  Basis  Set  Superposition  Effect,"  Int.  J.  Quantum  Chem.  23, 

847-854  (1983).  ” 

2.  P.C.  Hariharan,  W.S.  Koski,  Joyce  J.  Kaufman  and  Richard  S. 

Miller,  "Ab-Initio  MODPOT/VRDDO/MERGE  Calculations  on  Energetic 
Compounds  III.  Nitroexplosives :  Polyaminopolynitrobenzenes 

(Including  DATB,  TATB,  and  Tetryl),*"  Int.  J.  Quantum  Chem.  23,1493-1504  (1983). 

♦Paper  XVII  in  the  Series  "Molecular  Calculations  with  the  Non- 
Empirical  Ab-Initio  MODPOT/VRDDO/MERGE  Procedures." 

3.  W.  A.  Sokalski,  P.C.  Hariharan,  and  Joyce  J.  Kaufman,  "An  SCF 
Interaction  Energy  Decomposition  Study  of  12  Hydrogen  Bonded  Dimers," 

J.  Phys.  Chem.,  87,  2803  -  2810  (1983). 


B.  Accepted  for  Publication  and  In  Press 

1.  W.A.  Sokalski,  P.C.  Hariharan,  and  Joyce  J.  Kaufman,  "Guidelines 
for  Development  of  Basis  Sets  for  First  Order  Intermolecular  Energy 
Calculations,"  in  press,  J.  Comp.  Chem. 

2.  W.A.  Sokalski,  S.  Roszak,  P.C.  Hariharan,  Walter  S.  Koski, 

Joyce  J.  Kaufman,  A.H.  Lowrey  and  R.S.  Miller,  "Crystal  Structure 
Studies  Using  Ab-Initio  MODPOT/VRDDO  SCF  Energy  Calculations. 

I.  N2  and  CO2  Test  Cases.  II.  Nitromethane,"  an  invited  paper 
presented  at  the  Sanibel  International  Symposium  on  Molecular 
and  Solid  State  Theory  and  Computational  Quantum  Chemistry,  Palm 
Coast,  Florida,  March  1983.  In  Press,  Int.  J.  Quantum  Chem., 
Symposium  Issue. 

3.  P.C.  Hariharan,  Joyce  J.  Kaufman,  Alfred  H.  Lowrey  and  Richard 
S.  Miller,  "Ab-Initio  MODPOT/VRDDO/MERGE  Calculations  on  Energetic 
Compounds.  IV.  Nitrocubanes:  Mononitro-  to  Octanitro-  Quantum 
Chemical  Calculations  and  Electrostatic  Molecular  Potential  Contour 
Maps,"  an  invited  paper  presented  at  the  Sanibel  International 
Symposium  on  Molecular  and  Solid  State  Theory  and  Computational 
Quantum  Chemistry,  Palm  Coast,  Florida,  March  1983.  In  press,  Int. 

J.  Qunatum  Chem. 

4.  Joyce  J.  Kaufman,  "Ab-Initio  Potential  Energy  Functions  for 
Scattering  from  Energy  Partitioned  Ab-Initio  SCF  Energies  Plus 
Dispersion  -  Workshop  Note,"  presented  at  the  Sanibel  International 
Symposium  on  Molecular  Scattering  and  Quantum  Chemistry,  Palm  Coast, 
Florida,  March  1983.  In  press,  Int.  J.  Quantum  Chem.,  Symposium 
Issue. 
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X.  Papers  Published  on  This  Research  (cont.) 

B.  Papers  Accepted  for  Publication  and  In  Press  (cont.) 


5.  Cary  Chabalowski ,  P.C.  Hariharan,  Joyce  J.  Kaufman  and  Robert 
J.  Buenker,  "Ab-Initio  Multi  reference  Cl  Calculations  on  CH3N0,> 
Confirm  Earlier  Preliminary  GVB  and  MCSCF/CI  Results  that  HNO^ 
and  CH~N0?  Have  Multiconfiguration  Ground  as  Well  as  Electronically 
Excited  States  Even  at  Equilibrium  Geometry,  a  Symposium  Note," 
presented  at  the  Sanibel  International  Symposium  on  Molecular 
and  Solid  State  Theory  and  Computational  Quantum  Chemistry,  Palm 
Coast,  Florida,  March  1983.  In  Press,  Int.  J.  Quantum  Chem. , 
Symposium  Issue. 


XI.  Project  Personnel 


The  following  scientific  personnel  worked  full  or  part  time  with 
this  present  ONR  project. 

Joyce  J.  Kaufman,  Ph.D.,  Principal  Investigator 

P.  C.  Hariharan,  Ph.D.,  Associate  Research  Scientist 

Overall  responsibility  for  implementing  new  program  developments. 
Quantum  chemical  calculations  on  energetic  polymers. 

Quantum  chemical  calculations  on  the  large  nitroexplosives 
and  polynitropolyhedranes. 

Development  of  programs  for  graphical  display  of  three  dimensional 
isopctential  electrostatic  molecular  potential  contour  maps  for 
visual  display  and  for  computer  image  matching. 

Cary  Chabalowski,  Ph.D.,  Associate  Research  Scientist 

Program  implementation  of  large  scale  ab-initio  configuration 
interaction  program  based  on  multi  reference  determinants. 

Configuration  interaction  calculations  on  the  C-N02  decomposition 
pathway  of  nitromethane. 

Configuration  interaction  calculations  on  RDX. 

Terminated  September  15,  1983. 

W.  Andrezj  Sokalski,  Ph.D.,  Postdoctoral 

Further  development  and  implementation  of  program  for  partitioning 
of  the  SCF  wave  functions  and  energies  into  their  various  components 
and  addition  of  higher  order  terms. 

Carrying  out  SCF  intermolecular  calculations  for  additional  test 
molecules  for  CRYSTAL  calculations. 

Implementing  calculations  of  additional  terms  in  CRYSTAL  program 
itself. 

CRYSTAL  calculations  on  CH3NO2  and  RDX. 

Now  carrying  out  CRYSTAL  calculations  for  other  systems. 

Szczepan  Roszak,  Ph.D.,  Postdoctoral 

Implemented  program  for  many-body  and  coupled  cluster  multi  - 
determinant  repetitive  calculations  using  both  strictly  orthogonal ized 
localized  molecular  orbitals  and  conventional  canonical  delocalized 
molecular  orbitals. 

Coupled  cluster  test  calculations. 

Coupled  cluster  calculations  on  RDX. 

John  M.  Blaisdell,  Ph.D..  Postdoctoral 

Derived  and  implemented  ab-initio  crystal  orbital  program. 

Fitted  additional  potential  functions  for  molecules  for  CRYSTAL. 
Terminated  October  30,  1983. 
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Alexandre  Laforgue,  Ph.D.,  Visiting  Scientinst 

Derivation  of  new  approach  to  many-body  theory,  multideterminant 
reference  states  and  iterative  calculations. 

Terminated  December  1982. 
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